<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"><meta http-equiv="X-UA-Compatible" content="IE=edge,IE=9,chrome=1"><meta name="generator" content="MATLAB R2018b"><title>Untitled</title><style type="text/css">.rtcContent { padding: 30px; } .CodeBlock { background-color: #F7F7F7; margin: 10px 0 10px 0;}
.S0 { border-left: 0.833333313465118px solid rgb(233, 233, 233); border-right: 0.833333313465118px solid rgb(233, 233, 233); border-top: 0.833333313465118px solid rgb(233, 233, 233); border-bottom: 0.833333313465118px solid rgb(233, 233, 233); border-radius: 4px 4px 0px 0px; padding: 6px 45px 4px 13px; line-height: 17.234001159668px; min-height: 18px; white-space: nowrap; color: rgb(0, 0, 0); font-family: Menlo, Monaco, Consolas, 'Courier New', monospace; font-size: 14px;  }
.S1 { color: rgb(64, 64, 64); padding: 10px 0px 6px 17px; background: rgb(255, 255, 255) none repeat scroll 0% 0% / auto padding-box border-box; font-family: Menlo, Monaco, Consolas, 'Courier New', monospace; font-size: 14px; overflow-x: hidden; line-height: 17.234001159668px;  }
.embeddedOutputsErrorElement {min-height: 18px; max-height: 250px; overflow: auto;}
.embeddedOutputsErrorElement.inlineElement {}
.embeddedOutputsErrorElement.rightPaneElement {}
.embeddedOutputsWarningElement{min-height: 18px; max-height: 250px; overflow: auto;}
.embeddedOutputsWarningElement.inlineElement {}
.embeddedOutputsWarningElement.rightPaneElement {}
.diagnosticMessage-wrapper {font-family: Menlo, Monaco, Consolas, "Courier New", monospace; font-size: 12px;}
.diagnosticMessage-wrapper.diagnosticMessage-warningType {color: rgb(255,100,0);}
.diagnosticMessage-wrapper.diagnosticMessage-warningType a {color: rgb(255,100,0); text-decoration: underline;}
.diagnosticMessage-wrapper.diagnosticMessage-errorType {color: rgb(230,0,0);}
.diagnosticMessage-wrapper.diagnosticMessage-errorType a {color: rgb(230,0,0); text-decoration: underline;}
.diagnosticMessage-wrapper .diagnosticMessage-messagePart {white-space: pre-wrap;}
.diagnosticMessage-wrapper .diagnosticMessage-stackPart {white-space: pre;}
.embeddedOutputsTextElement,.embeddedOutputsVariableStringElement {white-space: pre; word-wrap: initial; min-height: 18px; max-height: 250px; overflow: auto;}
.textElement,.rtcDataTipElement .textElement {padding-top: 3px;}
.embeddedOutputsTextElement.inlineElement,.embeddedOutputsVariableStringElement.inlineElement {}
.inlineElement .textElement {}
.embeddedOutputsTextElement.rightPaneElement,.embeddedOutputsVariableStringElement.rightPaneElement {min-height: 16px;}
.rightPaneElement .textElement {padding-top: 2px; padding-left: 9px;}
.S2 { border-left: 0.833333313465118px solid rgb(233, 233, 233); border-right: 0.833333313465118px solid rgb(233, 233, 233); border-top: 0.833333313465118px solid rgb(233, 233, 233); border-bottom: 0px none rgb(0, 0, 0); border-radius: 0px; padding: 6px 45px 0px 13px; line-height: 17.234001159668px; min-height: 18px; white-space: nowrap; color: rgb(0, 0, 0); font-family: Menlo, Monaco, Consolas, 'Courier New', monospace; font-size: 14px;  }
.S3 { border-left: 0.833333313465118px solid rgb(233, 233, 233); border-right: 0.833333313465118px solid rgb(233, 233, 233); border-top: 0px none rgb(0, 0, 0); border-bottom: 0px none rgb(0, 0, 0); border-radius: 0px; padding: 0px 45px 0px 13px; line-height: 17.234001159668px; min-height: 18px; white-space: nowrap; color: rgb(0, 0, 0); font-family: Menlo, Monaco, Consolas, 'Courier New', monospace; font-size: 14px;  }
.S4 { border-left: 0.833333313465118px solid rgb(233, 233, 233); border-right: 0.833333313465118px solid rgb(233, 233, 233); border-top: 0px none rgb(0, 0, 0); border-bottom: 0.833333313465118px solid rgb(233, 233, 233); border-radius: 0px; padding: 0px 45px 4px 13px; line-height: 17.234001159668px; min-height: 18px; white-space: nowrap; color: rgb(0, 0, 0); font-family: Menlo, Monaco, Consolas, 'Courier New', monospace; font-size: 14px;  }
.S5 { border-left: 0.833333313465118px solid rgb(233, 233, 233); border-right: 0.833333313465118px solid rgb(233, 233, 233); border-top: 0.833333313465118px solid rgb(233, 233, 233); border-bottom: 0.833333313465118px solid rgb(233, 233, 233); border-radius: 0px; padding: 6px 45px 4px 13px; line-height: 17.234001159668px; min-height: 18px; white-space: nowrap; color: rgb(0, 0, 0); font-family: Menlo, Monaco, Consolas, 'Courier New', monospace; font-size: 14px;  }
.S6 { margin: 3px 10px 5px 4px; padding: 0px; line-height: 20px; min-height: 0px; white-space: pre-wrap; color: rgb(60, 60, 60); font-family: Helvetica, Arial, sans-serif; font-style: normal; font-size: 20px; font-weight: bold; text-align: left;  }
.S7 { margin: 20px 10px 5px 4px; padding: 0px; line-height: 20px; min-height: 0px; white-space: pre-wrap; color: rgb(60, 60, 60); font-family: Helvetica, Arial, sans-serif; font-style: normal; font-size: 20px; font-weight: bold; text-align: left;  }
.S8 { border-left: 0.833333313465118px solid rgb(233, 233, 233); border-right: 0.833333313465118px solid rgb(233, 233, 233); border-top: 0.833333313465118px solid rgb(233, 233, 233); border-bottom: 0px none rgb(0, 0, 0); border-radius: 4px 4px 0px 0px; padding: 6px 45px 0px 13px; line-height: 17.234001159668px; min-height: 18px; white-space: nowrap; color: rgb(0, 0, 0); font-family: Menlo, Monaco, Consolas, 'Courier New', monospace; font-size: 14px;  }
.S9 { border-left: 0.833333313465118px solid rgb(233, 233, 233); border-right: 0.833333313465118px solid rgb(233, 233, 233); border-top: 0px none rgb(0, 0, 0); border-bottom: 0.833333313465118px solid rgb(233, 233, 233); border-radius: 0px 0px 4px 4px; padding: 0px 45px 4px 13px; line-height: 17.234001159668px; min-height: 18px; white-space: nowrap; color: rgb(0, 0, 0); font-family: Menlo, Monaco, Consolas, 'Courier New', monospace; font-size: 14px;  }</style></head><body><div class = rtcContent><div class="CodeBlock"><div class="inlineWrapper outputs"><div  class = 'S0'><span style="white-space: pre;"><span>syms </span><span style="color: rgb(160, 32, 240);">clear</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsWarningElement" uid="A229D8D0" data-testid="output_0" data-width="559" data-height="30" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="diagnosticMessage-wrapper diagnosticMessage-warningType"><div class="diagnosticMessage-messagePart">Warning: SYMS X CLEAR and SYMS('X','CLEAR') will not accept CLEAR in a future release. SYMS X clears assumptions by default.</div><div class="diagnosticMessage-stackPart"></div></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'><span style="white-space: pre;"><span>clear</span></span></div></div><div class="inlineWrapper"><div  class = 'S3'></div></div><div class="inlineWrapper"><div  class = 'S3'><span style="white-space: pre;"><span>syms </span><span style="color: rgb(160, 32, 240);">z</span></span></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>N_1 = 1 - 3 * z^2 + 2 * z^3</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="948986AC" data-testid="output_1" data-width="559" data-height="29" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="symbolicElement"><span class="embeddedOutputsVariableElement">N_1 =&nbsp;</span><span class="MathEquation inlineSymbolicElement" style="font-size: 15px;"><span style="vertical-align: -5px;"><img src="" width="87" height="19"></span></span></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'><span style="white-space: pre;"><span>figure</span></span></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>fplot(N_1, [0, 1])</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsFigure" uid="B05AAB28" data-testid="output_2" style="width: 589px;"><div class="figureElement"><div class="figureContainingNode" style="width: 560px; max-width: 100%; display: inline-block;"><div class="GraphicsView" data-dojo-attach-point="graphicsViewNode,backgroundColorNode" id="uniqName_197_0" widgetid="uniqName_197_0" style="width: 100%; height: auto;"><img class="ImageView figureImage" data-dojo-attach-point="imageViewNode" draggable="false" ondragstart="return false;" id="uniqName_197_2" widgetid="uniqName_197_2" src="" style="width: 100%; height: auto;"></div></div></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>N_3 = 3 * z^2 - 2 * z^3</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="A760CBC0" data-testid="output_3" data-width="559" data-height="29" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="symbolicElement"><span class="embeddedOutputsVariableElement">N_3 =&nbsp;</span><span class="MathEquation inlineSymbolicElement" style="font-size: 15px;"><span style="vertical-align: -5px;"><img src="" width="62.5" height="19"></span></span></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'><span style="white-space: pre;"><span>figure</span></span></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>fplot(N_3, [0, 1])</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsFigure" uid="3C877CA5" data-testid="output_4" style="width: 589px;"><div class="figureElement" style="cursor: default;"><div class="figureContainingNode" style="width: 560px; max-width: 100%; display: inline-block;"><div class="GraphicsView" data-dojo-attach-point="graphicsViewNode,backgroundColorNode" id="uniqName_197_3" widgetid="uniqName_197_3" style="width: 100%; height: auto;"><img class="ImageView figureImage" data-dojo-attach-point="imageViewNode" draggable="false" ondragstart="return false;" id="uniqName_197_5" widgetid="uniqName_197_5" src="" style="width: 100%; height: auto;"></div></div></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'></div></div><div class="inlineWrapper"><div  class = 'S3'><span style="white-space: pre;"><span>syms </span><span style="color: rgb(160, 32, 240);">L</span></span></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>N_2 = L * z * (z - 1)^2</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="66519413" data-testid="output_5" data-width="559" data-height="29" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="symbolicElement"><span class="embeddedOutputsVariableElement">N_2 =&nbsp;</span><span class="MathEquation inlineSymbolicElement" style="font-size: 15px;"><span style="vertical-align: -5px;"><img src="" width="74" height="20"></span></span></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'><span style="white-space: pre;"><span>figure</span></span></div></div><div class="inlineWrapper"><div  class = 'S3'><span style="white-space: pre;"><span style="color: rgb(34, 139, 34);">% MuPAD code that caused issue 1: plot::Function2d(subs(N_2, L = 1), z = 0..1, L = -5..5)</span></span></div></div><div class="inlineWrapper"><div  class = 'S3'><span style="white-space: pre;"><span>aux1 = subs([0, 1], L, sym(5));</span></span></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>fplot(subs(subs(N_2, L, sym(1)), L, sym(5)), double(aux1))</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsFigure" uid="74C8EAB6" data-testid="output_6" style="width: 589px;"><div class="figureElement" style="cursor: default;"><div class="figureContainingNode" style="width: 560px; max-width: 100%; display: inline-block;"><div class="GraphicsView" data-dojo-attach-point="graphicsViewNode,backgroundColorNode" id="uniqName_197_6" widgetid="uniqName_197_6" style="width: 100%; height: auto;"><img class="ImageView figureImage" data-dojo-attach-point="imageViewNode" draggable="false" ondragstart="return false;" id="uniqName_197_8" widgetid="uniqName_197_8" src="" style="width: 100%; height: auto;"></div></div></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>N_4 = -(L * (z^2 - z^3))</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="97A3DBB6" data-testid="output_7" data-width="559" data-height="29" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="symbolicElement"><span class="embeddedOutputsVariableElement">N_4 =&nbsp;</span><span class="MathEquation inlineSymbolicElement" style="font-size: 15px;"><span style="vertical-align: -5px;"><img src="" width="79" height="20"></span></span></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'><span style="white-space: pre;"><span>figure</span></span></div></div><div class="inlineWrapper"><div  class = 'S3'><span style="white-space: pre;"><span style="color: rgb(34, 139, 34);">% MuPAD code that caused issue 2: plot::Function2d(subs(N_4, L = 1), z = 0..1, L = -5..5)</span></span></div></div><div class="inlineWrapper"><div  class = 'S3'><span style="white-space: pre;"><span>aux2 = subs([0, 1], L, sym(5));</span></span></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>fplot(subs(subs(N_4, L, sym(1)), L, sym(5)), double(aux2))</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsFigure" uid="D69A0DF5" data-testid="output_8" style="width: 589px;"><div class="figureElement" style="cursor: default;"><div class="figureContainingNode" style="width: 560px; max-width: 100%; display: inline-block;"><div class="GraphicsView" data-dojo-attach-point="graphicsViewNode,backgroundColorNode" id="uniqName_197_9" widgetid="uniqName_197_9" style="width: 100%; height: auto;"><img class="ImageView figureImage" data-dojo-attach-point="imageViewNode" draggable="false" ondragstart="return false;" id="uniqName_197_11" widgetid="uniqName_197_11" src="" style="width: 100%; height: auto;"></div></div></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'></div></div><div class="inlineWrapper"><div  class = 'S3'><span style="white-space: pre;"><span>syms </span><span style="color: rgb(160, 32, 240);">X</span></span></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>N_1 = subs(N_1, z, X/L)</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="E030C7A2" data-testid="output_9" data-width="559" data-height="64" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="symbolicElement"><div class="embeddedOutputsVariableElement">N_1 =&nbsp;</div><span class="MathEquation displaySymbolicElement" style="font-size: 15px;"><span style="vertical-align: -16px;"><img src="" width="96.5" height="36.5"></span></span></div></div></div></div><div class="inlineWrapper outputs"><div  class = 'S5'><span style="white-space: pre;"><span>N_2 = subs(N_2, z, X/L)</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="291AAEA3" data-testid="output_10" data-width="559" data-height="67" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="symbolicElement"><div class="embeddedOutputsVariableElement">N_2 =&nbsp;</div><span class="MathEquation displaySymbolicElement" style="font-size: 15px;"><span style="vertical-align: -15px;"><img src="" width="77" height="40"></span></span></div></div></div></div><div class="inlineWrapper outputs"><div  class = 'S5'><span style="white-space: pre;"><span>N_3 = subs(N_3, z, X/L)</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="D697CE07" data-testid="output_11" data-width="559" data-height="64" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="symbolicElement"><div class="embeddedOutputsVariableElement">N_3 =&nbsp;</div><span class="MathEquation displaySymbolicElement" style="font-size: 15px;"><span style="vertical-align: -16px;"><img src="" width="72.5" height="36.5"></span></span></div></div></div></div><div class="inlineWrapper outputs"><div  class = 'S5'><span style="white-space: pre;"><span>N_4 = subs(N_4, z, X/L)</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="5E748FBB" data-testid="output_12" data-width="559" data-height="66" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="symbolicElement"><div class="embeddedOutputsVariableElement">N_4 =&nbsp;</div><span class="MathEquation displaySymbolicElement" style="font-size: 15px;"><span style="vertical-align: -16px;"><img src="" width="99" height="38.5"></span></span></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>N = [N_1; N_2; N_3; N_4]</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="8F133084" data-testid="output_13" data-width="559" data-height="191" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="symbolicElement"><div class="embeddedOutputsVariableElement">N =&nbsp;</div><span class="MathEquation displaySymbolicElement" style="font-size: 15px;"><span style="vertical-align: -77px;"><img src="" width="125" height="163.5"></span></span></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span class="warning_squiggle_rte">B</span><span> = symhold(</span><span style="color: rgb(160, 32, 240);">'transpose(diff(diff(N, X), X))'</span><span>)</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="A99FDAA4" data-testid="output_14" data-width="559" data-height="67" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="symbolicElement"><div class="embeddedOutputsVariableElement">B =&nbsp;</div><span class="MathEquation displaySymbolicElement" style="font-size: 15px;"><span style="vertical-align: -15px;"><img src="" width="84" height="40"></span></span></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>B = transpose(diff(diff(N, X), X))</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="9DF6235B" data-testid="output_15" data-width="559" data-height="82" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="symbolicElement"><div class="embeddedOutputsVariableElement">B =&nbsp;</div><span class="MathEquation displaySymbolicElement" style="font-size: 15px;"><span style="vertical-align: -17px;"><img src="" width="380.5" height="54"></span></span></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'></div></div><div class="inlineWrapper"><div  class = 'S3'><span style="white-space: pre;"><span>syms </span><span style="color: rgb(160, 32, 240);">K</span></span></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>K </span><span class="warning_squiggle_rte warningHighlight">==</span><span> symhold(</span><span style="color: rgb(160, 32, 240);">'EI*int(transpose(B)*B, X, 0, L)'</span><span>)</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="B3B56395" data-testid="output_16" data-width="559" data-height="62" data-hashorizontaloverflow="false" style="width: 589px; max-height: 261px;"><div class="symbolicElement"><div class="embeddedOutputsVariableElement">ans =&nbsp;</div><span class="MathEquation displaySymbolicElement" style="font-size: 15px;"><span style="vertical-align: -12px;"><img src="" width="110.5" height="34"></span></span></div></div></div></div><div class="inlineWrapper"><div  class = 'S2'></div></div><div class="inlineWrapper"><div  class = 'S3'><span style="white-space: pre;"><span>syms </span><span style="color: rgb(160, 32, 240);">EI</span></span></div></div><div class="inlineWrapper outputs"><div  class = 'S4'><span style="white-space: pre;"><span>K = EI * int(transpose(B) * B, X, sym(0), L)</span></span></div><div  class = 'S1'><div class="inlineElement eoOutputWrapper embeddedOutputsSymbolicElement" uid="1C525C18" data-testid="output_17" data-width="559" data-height="280" data-hashorizontaloverflow="false" style="width: 589px; max-height: 291px;"><div class="symbolicElement"><div class="embeddedOutputsVariableElement">K =&nbsp;</div><span class="MathEquation displaySymbolicElement" style="font-size: 15px;"><span style="vertical-align: -123px;"><img src="" width="204.5" height="252.5"></span></span></div></div></div></div></div><h2  class = 'S6'><span>Local Functions</span></h2><h2  class = 'S7'><span>Helpers</span></h2><div class="CodeBlock"><div class="inlineWrapper"><div  class = 'S8'><span style="white-space: pre;"><span style="color: rgb(0, 0, 255);">function </span><span>h = symhold(exprstring)</span></span></div></div><div class="inlineWrapper"><div  class = 'S3'><span style="white-space: pre;"><span>h = evalin(symengine, [</span><span style="color: rgb(160, 32, 240);">'hold(' </span><span>exprstring </span><span style="color: rgb(160, 32, 240);">' )'</span><span>]);</span></span></div></div><div class="inlineWrapper"><div  class = 'S9'><span style="white-space: pre;"><span style="color: rgb(0, 0, 255);">end</span></span></div></div></div></div><br>
<!-- 
##### SOURCE BEGIN #####
syms clear
clear

syms z
N_1 = 1 - 3 * z^2 + 2 * z^3
figure
fplot(N_1, [0, 1])

N_3 = 3 * z^2 - 2 * z^3
figure
fplot(N_3, [0, 1])

syms L
N_2 = L * z * (z - 1)^2
figure
% MuPAD code that caused issue 1: plot::Function2d(subs(N_2, L = 1), z = 0..1, L = -5..5)
aux1 = subs([0, 1], L, sym(5));
fplot(subs(subs(N_2, L, sym(1)), L, sym(5)), double(aux1))

N_4 = -(L * (z^2 - z^3))
figure
% MuPAD code that caused issue 2: plot::Function2d(subs(N_4, L = 1), z = 0..1, L = -5..5)
aux2 = subs([0, 1], L, sym(5));
fplot(subs(subs(N_4, L, sym(1)), L, sym(5)), double(aux2))

syms X
N_1 = subs(N_1, z, X/L)
N_2 = subs(N_2, z, X/L)
N_3 = subs(N_3, z, X/L)
N_4 = subs(N_4, z, X/L)

N = [N_1; N_2; N_3; N_4]

B = symhold('transpose(diff(diff(N, X), X))')

B = transpose(diff(diff(N, X), X))

syms K
K == symhold('EI*int(transpose(B)*B, X, 0, L)')

syms EI
K = EI * int(transpose(B) * B, X, sym(0), L)
%% Local Functions
%% Helpers
%%
function h = symhold(exprstring)
h = evalin(symengine, ['hold(' exprstring ' )']);
end
##### SOURCE END #####
--></body></html>